Fission Barriers of Compound Superheavy Nuclei 
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The dependence of fission barriers on the excitation energy of the compound nucleus impacts 
the survival probability of superheavy nuclei synthesized in heavy-ion fusion reactions. In this 
work, we investigate the isentropic fission barriers by means of the self-consistent nuclear density 
functional theory. The relationship between isothermal and isentropic descriptions is demonstrated. 



Calculations have been carried out for Fm, 



^Ds, 



'112, 



114, and ^^^124. For nuclei around 



112 produced in "cold fusion" reactions, we predict a more rapid decrease of fission barriers with 
excitation energy as compared to the nuclei around ^^^114 synthesized in "hot fusion" experiments. 
This is explained in terms of the difference between the ground-state and saddle-point temperatures. 
The effect of the particle gas is found to be negligible in the range of temperatures studied. 



PACS numbers: 24.75.-(-i, 21.60.Jz, 27.90.-fb, 24.10. Pa 

What are the heaviest nuclei that can exist? To an- 
swer this question, nuclear physicists explore superheavy 
systems at the limit of mass and charge. During recent 
years, the field has witnessed remarkable progress tl„i3j3 
in the production and identification of new elements. The 
major experimental challenge is to find optimal beam- 
target combinations and kinematic conditions that would 
lead to the formation, at reasonable rates, of the species 
of interest. One of the key problems is the survival prob- 
ability of a superheavy nucleus synthesized in a heavy-ion 
fusion reaction that depends on a competition between 
fission and particle evaporation 

The dependence of the fission barrier on the excitation 
energy is among the key factors determining the pro- 
duction of a superheavy nucleus. Since shell effects, es- 
sential for the mere existence of superheavy nuclei, are 
quenched at high temperatures (see, e.g., Refs. 0, §]), 
it is expected that the fission barriers in superheavy CN 
should quickly decrease with energy. In the analysis of 
experimental data, this is usually accounted for by a phe- 
nomenological damping factor 0, 0, H, Q (cf. discussion 
in Ref. [i]). 

Microscopically, shell effects in superconducting 
heated nuclei can be self-consistently treated by the 
Finite- TemperatureHartree-Fock-Bogoliubov (FT-HFB) 



method 

SHI [ll|, For superheavy nuclei, although 
there have been extensive self-consistent studies of zero- 
temperature fission pathways, (see, e.g., [l^, [3) fisj)- 
studies of CN fission have been virtually nonexistent. 
Moreover, in the majority of studies, CN fission has been 
treated as an isothermal process in terms of free energy, 
F=E — TS at a fixed temperature T. The assumption 
of r=const is certainly wrong: the fissioning CN is not 
connected to an external thermal reservoir. Physically, 
a more appropriate picture of fission is that governed by 



the isentropic process [16|. The necessary condition for 
the isentropic scenario is that the collective motion is 
adiabatic, i.e., no heat energy can be delivered to nor 
extracted from the system. 

The aim of this study is to investigate self-consistent 
isentropic fission pathways in superheavy CN, with a fo- 
cus on energy dependence of fission barriers. To this 
end, we selected several nuclei of current experimental 
interest: (i) ^^^Ds and ^"^^112* that have been synthe- 
sized in the "cold-fusion" reaction using a ^°*Pb tar- 



208j 

get at excitation energies E* of ~10-12 MeV 0; (n) 
the nucleus '^^'^114* produced in the "hot-fusion" reac- 
tion 48Ca-|-244pu at £;*-36-40MeV [i3l; (iii) the nucleus 
312124* at £;*-80MeV studied by means of '''^Ge+^^^V 
reaction and crystal blocking [18|; and (iv) the ^^'^Fm 
that is expected to fission symmetrically into two doubly- 
magic Sn nuclei. The FT-HFB calculations are carried 
out using the two Skyrme-HFB codes: the rece ntly de- 
veloped axial coordinate-space solver HFB-AX [19| and 
a symmetry- unrestricted solver HFODD The de- 

scription of thermal properties involves significant contri- 
butions from high-lying single-particle states which giv e 
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rise to the particle gas as the temperature increases 
requiring a very large configuration space to guarantee 
convergence j22!|. In this respect, HFB-AX is an excel- 
lent tool as it allows calculations in very large deformed 
boxes. The finite-temperature formalism has been imple- 
mented in HFB-AX and HFODD in the usual way [11 1 
by introducing the thermal-averaged particle and pairing 
densities through the Fermi distribution function. 

In the particle-hole channel, we use the SkM* en- 
ergy density functional f23| that has been optimized at 
large deformations; hence it is often used for fission bar- 
rier predictions. In the pairing channel, we adopted 
the density-dependent 6 interaction in the mixed variant 
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[2i|. The pairing strengths, Vp=-332.5 MeVfm-^ (pro- 
tons) and 14,=— 268.9 MeVfm"'^ (neutrons) have been 
fitted to reproduce the pairing gaps in ■^'"'^Fm. The de- 
tails of HFB-AX calculations follow Ref. [ig. We used 
M—13 order B-splines, and the maximum mesh size 
h—0.6 fm. The cylindrical box employed depends on the 
total quadrupole moment of the system, Q20 ■ That is, for 
Q20<30b we used a square box of Rp=Rz 



=20.4 fm; for 
21.6 fm; and 
The 

calculations with HFODD were carried out in a space 
of the lowest 1161 stretched oscillator states originating 
from the 31 principal oscillator shells. 



30<Q2o<80b we took i?p=19.2fm and i?^ 
for Q20>80b we took i?p=18fm and i?^=22.8fm 
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FIG. 1: (Color online) Calculated isothermal [F at constant 
T) and isentropic [E at constant 5") axial symmetric fission 
pathways for ^^^114 as a function of the total quadrupole 
moment Q20- The isothermal calculations were carried out 
at fcT=0.5 and 1.5 MeV. In the isentropic description, 5* was 
fixed at the free energy minimum, i.e., S'/fc=18.5 (69.2) at 
/fcr=0.5 (1.5) MeV. 

As pointed out in Ref. ^li*!, isothermal and isentropic 
descriptions can be related by making use of the thermo- 

dynamical identity (^^^j — {^'S§^ which is remi- 
niscent of the well-known relation (^)s = (f^)^' 
deed, Q20 enters the variations dF and dE via the term 
—<l2QdQ20i where 920 is the Lagrange multiplier corre- 
sponding to the constraint on the quadrupole moment. 
Figure [T] displays isothermal and isentropic axial fission 
pathways of ^^^114 as a function of Q2o- In the isentropic 
description, S was fixed at the value S = S{T) corre- 
sponding to the free energy minimum at temperature T. 
This was done by performing constrained FT-HFB cal- 
culations for a number of temperatures and inverting the 
relation S = S(T) numerically by using interpolation. 
It is seen that the isothermal and isentropic curves are 
very close. It is worth noting that in the macroscopic- 
microscopic calculations of Ref. [l^ the isentropic bar- 
riers are predicted to be significantly higher than the 
isothermal ones. The reason for this is the violation of 
self-consistency in the macroscopic-macroscopic theory. 
Figure [1] shows that the variational principle behind FT- 



HFB guarantees practical equivalence of isothermal and 
isentropic pictures. The remaining small discrepancy is 
due to the numerical interpolation error caused by ex- 
traction of {E, S) values from the original {F, T) mesh. 
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FIG. 2: (Color online) The temperature along the isentropic 
symmetric fission pathways of ^^*112 and ^^^114. The ground 
state (g.s.) and saddle point configurations are marked by 
stars. The g.s. temperature was assumed to be T^.s. ~750 keV 
and S = S{Tg.s.). The temperature at the saddle point is 
considerably lower than Tg,a.. 

The behavior of temperature T = T{S) is shown 
in Fig. [2] along the fission pathways of ^^^112 and 
292^24. The entropy corresponds to the g.s. value at 
Tg,s.^750keV. It is seen that T changes as a function 
of Q20- In particular, the barrier temperature is signif- 
icantly lower than Tg^.- In the following, we shall stick 
to the isentropic description, i.e., the temperature will be 
related to the g.s. excitation energy. 
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FIG. 3: (Color online) Symmetric isentropic fission pathways 
of ^^*Fm at the values of S corresponding to kTg, s.=0, 0.5, 
1.0, and 1.5 MeV The energy has been normalized to zero at 
the ground-state minimum. The effect of triaxial degrees of 
freedom on the first barrier is marked by dashed lines. 

Figure [3] displays the energy curves of ^^''Fm as func- 
tions of Q20 at different entropies corresponding to dif- 



ferent values of Tg.s.. At S{kTg, 



0.5) , the barrier 
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increases by about 0.6 MeV as compared to S=0, due to 
the reduction of pairing correlations @, [l^, |Tl|, . (In 
our calculations, pairing energies are unimportant above 
fcT=0.7MeV.) At higher excitation energies, the barrier 
is gradually reduced to 0.9 MeV at S{kTg,s. — 1-5), due 
to the thermal quenching of shell effects. In order to 
estimate the reduction of the fission barrier due to triax- 
iality expected in the Fm isotopes l^, 2^ , we performed 
symmetry- unrestricted calculations with HFODD. The 
result is shown in Fig. [3] by a dashed line. At low exci- 
tation energies, triaxiality reduces the fission barrier by 
^^3-4 MeV, but the triaxial shell effect is washed out with 
increasing entropy and becomes negligible at the largest 
excitations considered. For the systematic calculations 
of triaxial and reflection-asymmetric deformations along 
the isentropic fission pathways of superheavy nuclei, we 
refer the reader to Ref. [26]. 
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FIG. 4: (Color online) The height of the inner fission bar- 
rier in ^'^^Fm, ^''^Ds, ^''^112, ^^^114, and ^^^124 as a function 
of excitation energy E* . The effect of triaxiality on fission 
barrier has been included. The experimental values of E* 
corresponding to CN formed in the reactions indicated are 
marked by arrows. 

To study the excitation energy dependence of fission 
barriers in more detail, in Fig. 2] we plot the height of 
the inner axial fission barrier of ^^"^Fm, ^^^Ds, ^^^112, 
^^^114, and ^^^124 as a function of the excitation en- 
ergy E*=E{S) ~ E(S = 0). Above fcTg.^ =0.5MeV, 
fission barriers Eb are damped exponentially with E*: 



to shell effects. As seen in Fig. [2l in the isentropic 
picture, the saddle point temperature Tg is lower than 
Egg , i.e., AT=Tgs — Tb>0 (in a nice analogy to 
an adiabatically expanding gas). Due to shell ef- 
fects, AT(278ll2)>Ar(292ii4) and the thermodynam- 
ical identity {^)q^^ ~ kT implies a larger in ^^*112, 
thus explaining the pattern seen in Fig. [3) 
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FIG. 5: (Color online) Proton (thick lines) and neutron (thin 
lines) spherical density distributions in ^^^114 calculated at 
different temperatures (in MeV). The external gas contribu- 
tions are marked by dotted lines. The number of nucleons 
in the gas is indicated by numbers. The size of the box in 
HFB-AX was 20.4 fm. 



As discussed in Ref. [27|, the FT-HFB solution in a 
confined box (or in a finite localized basis) corresponds 
to a nucleus located at the center surrounded by the ex- 
ternal gas. The gas produces the pressure necessary to 
obtain an equilibrium with the particle-decaying hot nu- 
cleus; the corresponding particle-decay lifetime is in fact 
inversely proportional to the density of the external gas 
2l| . Figure [5] illustrates the effect of the gas for spheri- 
cal configuration in '^^'^114 at different temperatures. In 
order to separategas contributions, we applied the pro- 



soce 



cedure of Refs. [21|, |28|. It is seen that with increasing 
temperature, the gas gradually appears. As expected, 
the neutron gas is uniformly distributed within the vol- 
ume of the box while the proton gas appears outside the 
nuclear surface due to the effect of the Coulomb bar- 
rier. The contribution from the gas to the kinetic en- 
ergy is very small: it is about 0.2 (0.9) MeV at /cT=I.25 
(1.5) MeV. The number of gas neutrons is very small; 



The value of the damping parameter it increases from A^gas=0.01 at fcT^lMeV to 0.21 at 



is not known well (see discussion in Ref. [9|) but is 
usually taken in the range of 8-20 MeV • According 
to our calculations, in the CN ^^^Ds and ^^^112 synthe- 
sized in cold fusion reaction 7^^ = 17.2 and 10.8 MeV, re- 
spectively, while it is 30 MeV in 292114 and 20.2 MeV in 
312124. 

The appreciable change in with N and Z, e.g., 
when going from ^"^^112 to ^^^114, can be traced back 



fcT=1.5MeV, cf. Fig.O Consequently, the gas contribu- 
tion to the deformation energy is practically negligible in 
the range of temperatures considered. At very high tem- 
peratures, however, the gas is expected to significantly 
change the Coulomb energy, the total entropy, and the 
chemical potentials, and its contribution should be prop- 
erly removed 21 1 . 

In conclusion, we performed self-consistent calcula- 
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tions of isentropic fission barriers of compound su- 
perheavy nuclei based on a coordinate-space FT-HFB 
method. We first demonstrated the relationship between 
the isothermal and isentropic description of fission and 
emphasized the role of self-consistency. The conclusion, 
important for practical applications, is that the surfaces 
of F(T = To) and E{S = So) in the space of collec- 
tive coordinates are identical for 5*0 = S{To) if the self- 
consistency condition is met. It is to be noted, however, 
that this formal connection does not indicate that the 
isothermal and isentropic pictures of fission are similar. 
The isothermal approach to fission is clearly incorrect as 
the nucleus is not connected to a thermal bath, i.e., the 
temperatures of the g.s. configuration and the barrier 
obviously differ. On the other hand, if the fission process 
is adiabatic, the isentropic description should be closer 
to reality. 

Secondly, we demonstrate that the dependence of 
isentropic fission barriers on excitation energy changes 
rapidly with particle number, pointing at the importance 
of shell effects even at large excitation energies charac- 
teristic of compound nuclei. For instance, fission barriers 
for ^^^Ds and ^^^112, produced in a cold fusion reaction, 
and ^^^114, synthesized in a hot fusion reaction, are pre- 
dicted to exhibit markedly different behavior. For the 
CN ^^^124, we calculate no isentropic fission barrier at 
£;*-80MeV. 

Finally, we show that the external particle gas has no 
effect on the fission barriers up to at least A;T=1.5MeV. 
The barrier damping parameters extracted from our FT- 
HFB calculations, as well as the neutron decay rates ex- 
tracted from the magnitude of the neutron gas compo- 
nent [2l| can be used to provide reliable theoretical esti- 
mates of CN survival probability. Work along these lines 
is in progress. 
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